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Abstract 

We give an independent proof of the Krasikov-Litsyn bound d/n < (1 — 5 _1//4 )/2 on doubly-even self-dual 
binary codes. The technique used (a refinement of the Mallows-Odlyzko-Sloane approach) extends easily 
to other families of self-dual codes, modular lattices, and quantum codes; in particular, we show that the 
Krasikov-Litsyn bound applies to singly-even binary codes, and obtain an analogous bound for unimodular 
lattices. We also show that in each case, our bound differs from the true optimum by an amount growing 
faster than 0(y / n). 

1 Introduction 

In |l(J, Q, Mallows, Odlyzko, and Sloane proved the following result: 

Theorem. Let C be a doubly-even binary self-dual code of length n and minimum distance d. Then d < 
4[n/24] + 4. Morever, for any constant b, one has d < n/6 — b for sufficiently large n. 

as well as analogous results for ternary codes and even unimodular lattices. The first claim has since been 
extended to singly-even binary self-dual codes jL5l , and analogous results have been obtained for even strongly 
modular lattices Jl2| , |lg| ] and odd strongly modular lattices [[19), including the odd unimodular case. 

Regarding the asymptotic claim, essentially the only improvement is the bound of Krasikov and Litsyn [Q : 

Theorem. Let Ci be a family of doubly-even binary self-dual codes of length tending to infinity. Then 

.. d{Cj) ^ 1 - 5-v 4 

hmsup — — < . (1.1) 

However, it is unclear to what extent their argument extends to the other cases of interest (especially the lattice 
cases). 

In the present paper, we give a new technique for deriving bounds on self-dual codes and modular lattices. 
In the doubly-even binary case, our bound is precisely the Krasikov-Litsyn bound; the difference is that our 
technique easily generalizes. The basic idea (following ||) is to use invariant theory to construct linear relations 
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that must be satisfied by the weight enumerator of a self-dual code. In the doubly-even binary case, the simplest 
such relation gives d < 4[n/24] + 4; to obtain their asymptotic improvement, Mallows, Odlyzko, and Sloane 
also take into account the second-simplest relation. Somewhat surprisingly, our bounds also only use these two 
relations; equivalently, we only consider the coefficients of the weight enumerator up to weight 4[rt/24] + 8. 

To compute the coefficients of these relations, we use the Biirmann-Lagrange theorem to express them as 
coefficients of certain power series, which we asymptotically analyze via Cauchy's integral and the saddle-point 
method. Under suitable conditions on the power series, we can then show that their coefficients in certain 
ranges are asymptotically uniformly positive, and thus give a contradiction unless the minimum distance bound 
holds. This necessitates a certain amount of analysis, which we deal with in Section |^. Then, in Section |3[ 



we give our main theorem, Theorem 3.6. This is stated in some generality (regarding the minimum valuations 



of certain families of power series with nonnegative coefficients), so as to include most of our applications as 



special cases. We state these special cases in Section m see Theorem 4.1 for (most) codes, Theorem 4.2 for 



even modular lattices, and Theorems 4.3 and |4.5| for codes over Z4. The remaining applications not directly 



dealt with by Theorem 3.6 are considered in Section ^ see Theorem [5.1| for (singly-even) self-dual binary codes, 



Theorem 5.2 for odd modular lattices, and Theorems |5.5| and 5.E for quantum codes. 

As we remarked, our improvements on the main Mallows-Odlyzko-Sloane bounds are obtained by considering 
the first two relations coming from invariant theory, rather than just the first. In Section ^, we consider the 
possibility of improving the bounds by using the first k relations. In fact, it turns out that, despite the 
significant improvement between k = 1 and k — 2, increasing k beyond 2 does not give a better bound on 



limsupd/n. We do, however, obtain a slight lower-order improvement (Theorem 3.1); increasing k gives an 
0(n -1 / 2 ) improvement on the bound on d/n. In particular, each of our bounds differs from the true optimal 
minimum distance (norm) by an amount growing faster than 0(^fn). 

Acknowledgements. The author would like to thank H. Landau, A. M. Odlyzko, and N. J. A. Sloane 
for helpful discussions regarding Section ||, especially Lemma |2.3[ as well as I. Duursma for pointing out that 
Krasikov and Litsyn had improved their earlier bound to the one stated above. 

2 Lemmas 

We will use the notation [P}f(t) to refer to the coefficient of V in the (formal) Laurent series fit). We also use 
the notation 

fhg (2.1) 

to say that [P](f(t) - g(t)) > for all j. 

Lemma 2.1. Let f(t) be a Laurent series convergent on an annulus < r\ < \t\ < r2 < 00. If f y 0, then we 
have the bound 

\f(t)\ < f(\t\) (2.2) 

valid on the annulus. If \t 3 ]f(t) and are both nonzero for some j , then equality can hold only when 

t = \t\. 
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Proof. We have: 

i/wih£/^'i<£/;I*i j ' = /(i*d< ( 2 - 3 ) 

3 3 

with equality only when there exists a such that V — a\t\ 3 for all j with fj nonzero. Dividing two consecutive 
such equations, we obtain t = \t\ as required. □ 

The Hadamard three-circles theorem then immediately implies that log(/(e s )) is strictly convex, for n < 
e s < r 2 . In fact, we have the slightly stronger statement: 

Lemma 2.2. let f(t) be a Laurent series convergent on an annulus < r\ < \t\ <r 2 <oo ; and not proportional 
to z n for any n. If f h 0, then for r\ < e s < r 2 , 

^log(/( e s ))>0. (2.4) 

Proof. Setting r = e s , we have: 

* w( ,,,-(^)W>,-^-(^). 

Now, f(r) is positive on n < r < r 2 , so we may freely multiply by /(r) 2 ; we thus need to show that 

d \ 2 ., . ( & 



f(r)(r-j f {r )-{r-f(r)) > 0. (2.6) 
Now, the left-hand-side has a Laurent series convergent in the original annulus, namely 

(£ fc2 / fc * fc ) - (£i/^ j )(£ fc /^ fe ) -EE fc 2 w +fe - EE^/^' +fc ( 2 - ? ) 

j k 3 k j k j k 



(2.8) 

Since this y 0, and has at least one nonzero coefficient, the desired inequality follows. □ 



j k 



Lemma 2.3. Let F,G be real power series both convergent in the circle \t\ < r . Suppose furthermore that 
G y 0, with both G(0) and G'(0) positive. Then for any compact subset I C [0, r ) on which F is positive, and 
for all sufficiently large n (depending on I): 

[t m ]F(t)G(t) n > (2.9) 

whenever m/n e S(I), with S(t) := tG'(t)/G(t). 
Proof. It will suffice to show that 

[t s(r ^ n }F{t)G(t) n = (F(r) + o(l))[i s W n ]G(i)™, (2.10) 
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with error uniform on any interval r € [0, a] with a < r , since then on /, F(r) is bounded away from 0, while 
the error converges uniformly to 0. We split into two cases: < r < n~ 2 / 3 and n -2 / 3 < r < a. 
In the first region, we claim that for all sufficiently large n, and for < r < n~ 2/<3 , 



\ t s(r)n ]F{t)G{tr = {F{Q) + Q(n -i/s ))G(())n l (S(r)n)\ 



with error uniform in r. By Cauchy's residue theorem, 



G(0)- n [t s ^ n ]F(t)G(tr = ±-.t t~ s ^ n F(t) f^Y 

27T« J\t\=S(r)/S>(0) V G (°)/ 



Now, \t\ = 0(n 2 / 3 ), so we have the uniform estimates 



and thus 



Now, 



so it remains to show that 



/|t|=S(r)/S'(0) 

is bounded. But, setting m = S(r)n and rescaling t, this is 



'1*1 

by the known asymptotics of Bessel functions. 



f e mKe {t) dt = m\I (m) = 
m m 2m y, t , =1 t m m w ; 



Since 



\og\G(re l6 )\ -logG(r) 
7(P 



(2.11) 



(2.12) 



log G(t) - logG(O) = S'(0)t + 0(n~ 4/3 ), (2.13) 
F(t) = F(0) + 0(n- 2/3 ), (2.14) 



FW ^G|y) n = eS ' (0) " t( ^ (0) + o( ^ 1/3)) - (2 - 15) 



— f t- s ^ n F(0)e s '^ nt dt/t = F(0) {nS '} 0) ^? n , (2-16) 

2tt* 7| t | = s ( r)/s'(o) (S(r)n)! 



(5(r)n)!(5'(0)n)^«"-L f | exp(S'(0)ni)|y (2.17) 



(2.18) 



We now consider the case n 2 / 3 < r < a. Here, we claim 

[t s ^ r)n }F{t)G{t) n = (2TTnrS'{r))- 1 / 2 r- s(r ^G(r) n (F(r) + O^ 1 / 60 )), (2.19) 
again with uniform error. Again, Cauchy's integral gives 

r s ^ n G{r)- n [t s( - r)n ]F{t)G(t) n = 7^ f e~ iS ^ ne F(re ie ) (^Jp) d0. (2.20) 



(2.21) 
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is continuous and negative for r £ [0, 6], 9 € [—ir, n], there exists a positive constant G such that 

log \G(re ie ) | - log G(r) < -Cr6 2 (2.22) 

in that region. 

In particular, when \9\ > n - 9 / 2 o r -i/2^ we g n( j nr g2 > n 1 / 1 ^^ an d thus 



G{re w ) 



G(r) 



< e- Cnl/1 °. (2.23) 



Since \F(re l9 )\ is bounded, we find that the contribution of this region to the integral is negligible. 
Now, consider the region |6*| < n _9 / 20 r -1 / 2 . Since 

log G(re ie ) - log G(r) - iS{r)6 + rS"(r)f 
is continuous for r € [0, 6], € [— 7r, 7r], we find the uniform estimate 
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Similarly, 



and thus 



Since 



1 / e — s 'W e2 / 2 ^ = (2^nr5'(r))- 1 / 2 + O( e -" 1/10s 'M/ 2 ), 

271" J|fl|<ri-9/2i)r- 1 /2 



(2.24) 



logG(re ie ) - logG(r) - iS(r)6 + rS '^° = 0(r6 3 ). (2.25) 
or upon exponentiation, 

e- iS{r)ne G{r)- n G{re i6 ) n = e^ rS ' W92/2 (l + 0(nr9 3 )) = e ~^s'(r)e 2 /2^ + Q^-i/eo^ q.2G) 



F{re w ) = F(r) + 0(n- 7 / 60 ), (2.27) 



F{re l9 )e- lS{r)ne G{r)- n G{re ie ) n = e- nrS ' {r ^/ 2 {F{r) + O^ 1 / 60 )). (2.28) 



(2.29) 



271" ,/|0|<7i-9/2O r -l/2 

the claim follows. □ 

Remark. Note that we only used the fact G h through the conclusions of the previous two lemmas. 
Away from 0, we can give much stronger estimates (which will be used in the final section): 

Lemma 2.4. Fix radii < r\ < r 2 < oo, and let F(t) and G(t) be Laurent series convergent on a neighborhood 
of the annulus r\ < \t\ < r 2 ; suppose further that G(t) y and has two consecutive nonzero coefficients. Let 
S(r) := rG'(r)/G(r), and for each n < r < r 2 , define a power series j(r,x) by 7(r, 0) = 0, 7 x (r, 0) > and 

logG(re^( r ^) = log G{r) + iS(r)-f(r, x) - rS'{r) — . (2.30) 
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Then for all integers k > and for r\ <r<r 2 , we have the asymptotic estimate 

^2^rS'{r)nr s ^ n G{rr n [t s ^ n ]F{t)G{t) n = £ ( r5 '0>)~ J (-f^ lx (r, x)F(re i ^)) + 0(n- fe ), 

o<j<fe-i z J - ^ ax ' ^=o 

(2.31) 

with error uniform in r. 
Proof. Note that 

v 2 

logG(re iy ) = logG(r) + iS(r)y - rS"(r)y + 0(y 3 ); (2.32) 

since rS'(r) > 0, we conclude that 7(7-, x) converges for \x\ < x for some x > independent of r, and satisfies 
7x(r,0) = 1. 

By Laurent's theorem, 

v/27rrS'(r)nr s M"G(r)-"[i s M"]F(t)G(f)" = \[^^ f F{re w ) (^—-^^-1^ d6. (2.33) 

Now, as before, we can restrict the integral to any uniform neighborhood of 0, with exponentially small error. 
In particular, we may restrict to a neighborhood \6\ < 9 affording the change of variable 9 — j(r,x). The 
integral thus becomes 

rS'(r)n ' "' 
2t 

Now, we have the uniform estimate 



■ f ° lx {r,x)F{re l ^ x) )e- nrS,( - r)x2/2 dx. (2.34) 

J —Xn 



Jx (r,x)F(re i ^) = £ -7^(^—j7x{r, x)F(re ij( - r,x ^)) x= o + 0(x 2k ); (2.35) 

0<j<2fc-l J ' x 

since 

y/rS'(r)n2n f" x 2k e~ nrS ' ^ ^ dx < ^rS'(r)n2n f°° x 2k e~ nrS ' ^ ^ dx = (rS 1 (r)n)- k ^ , (2.36) 

J-x J-00 2 K. 

the contribution of the error term is as required. Once we remove this term, the integral can be extended to 
00, again giving uniform exponentially small error. Evaluating the resulting Gaussian integral gives the desired 
result. □ 

For our purposes, we will need a version of this valid in the neighborhood of a zero of F. Define polynomials 
hk(x) for integers k > via the generating function 

£M*4 = e * x " t2/2 ; ( 2 -37) 

fe>0 



in particular, hk(x) is a monic polynomial of degree k, and 

'¥j\(k-2j)\ 



0<j<[k/2] 



]T [% k ~ j hM- (2-39) 



0<j<k 
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(Thus hk(x) are just rescaled Hermite polynomials.) In terms of hj(0), the above estimate becomes (replacing 
k by fe/2): 



V2irrS'(r)nr s ( r)n G(r)- n [t s{r)n ]F(t)G(t) n = 

^ (rS'{r)n)-^hM (* {r , x)F{n **T*)^ + O { n -0+W). (2.40) 
^a^u 3- \dxi / x =o 



o<j<fc 



Corollary 2.5. M^ii/i hypotheses as in Lemma 2.4, suppose F has a zero of order k at the point r £ [^i,^] 
Then we have the uniform asymptotic estimate 



^^.^ ^((^v^ ^^.,,,.^,,, ( , 41) 

mZid /or ri < r < r2 . 

Proof. Since the main term of the estimate has order O(\r — ro\ k ), we can tolerate a multiplicative error of order 
1 + 0(n -1 / 2 ). In particular, we may replace [t s ^ n ]G(t) n by its first-order estimate. We thus need to estimate 



(2.40) 



Now, 



7 x (r ) x)F(re i ^) = J2((r^) l F)(r) lx (r,xh(r,x) 1 ^ (2.42) 



0<( 



Thus we obtain 



^ ( ! ^rW^) !i ^ _t {{r ^y F)(r) f ( _ly+i y{ r -x y+A . (2 , 5) 



Now, for < I < k, 

d 



and for I < j, we have: 



(r^)'F(r) = O((r-r ) fe -'), (2.46) 



(±y+^(r, X ) l+1 ) =()il). (2.471 



In particular, the j, i term gives a contribution of order 0(\r — ro\ k ~ l n~i/ 2 ). For fixed /, the contributions get 
smaller as j increases. We thus find that the terms with I < j < k are dominated by the terms with j = I + 1, 
of order 

0(\r ~ r \ k - l n- {l+1 ^ 2 ) = 0(max(|r - r |, n" 1 / 2 )^ 1 ). (2.48) 
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It remains to consider the terms with j — I, that is, 

E ^'^''^h i^jFm (2.49) 
o<;<fc 

If we replace r-^ by f*o^, F(r) by F^(ro)(r — ro) k /k\ and rS'(r) by roS"(ro), the resulting error is again 
0(max(|r — r |, n _1 / 2 ) fc+1 ). We thus obtain 



k\ o ^ fc V {S<{ro)n/rv) l /*\l) K u; k\ (S'(r )n/r ) fe / 2 

as required. □ 
We conclude with one more analytical lemma, used in Section [| below. 

Lemma 2.6. Let F(t) and G(t) be real power series convergent on the circle \t\ < r. Suppose G >z 0, with 
G(0), G'(0) both positive. Let L be a compact subset of [0,r) on which F is positive. Then for all sufficiently 
large I, the holomorphic functions 

F(a)G {a y + F(-a)G(-ay and ^W^Z^h^ (2 . 51) 

a 

are positive on I . 

Proof. We first note that if L does not contain 0, then \ G(—a)/G(a) \ is bounded below 1 on I, while \F(—a)/F(a)\ 
is bounded; the result follows immediately. If / is the single point {0}, then F(0) > by assumption, so F(x) > 
for sufficiently small positive x] we may thus enlarge I while maintaining the hypotheses. We may thus take / 
of the form [0, e] with e > 0; moreover, if the theorem is true for [0, e'] with < e' < e, it is true for [0, e]. 

Since F(0) and G(0) are both positive, we may choose e so that F(a) and G(a) are positive on [— e, e]. In 
particular, the first function is thus positive on [0, e] for all I; it remains to consider the second function. 

Now, clearly (F(a) — F(—a))/a is analytic on [— e, e], while 

logG(a)-logG(-a) 
a 

is real analytic and positive on [— e, e]. In particular, there exist positive constants C\ and C2 such that 



(2.52) 



\F(a) - F(-a)\ < C x a (2.53) 
C 2 a < log G(a) - log G(-a). (2.54) 



We thus find: 



F( W -F(-,) G (- fl )' = G{ _ a yF(a)-F { -a) + ^^^Oja)-^ - ! ^ 

> G(-aY ( F ^)-H-a) + logG( )-logG(- a) \ 

\ a a J 

> G(-a) l (F{a)C 2 l - d), (2.57) 
which is clearly positive for sufficiently large I. □ 
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3 The main theorem 

Given a power series (or left-finite Laurent series) p, we define the valuation v{p) of p to be the exponent of the 
first nonzero term in p. 

The general scenario we consider is as follows. We are given power series f(t), g(t), and h(t), and asked to 
prove a statement of the following form. If 

A(t) = h(t) Yl cifar-'gity (3.1) 

0<i<m 

is such that A >z 0, then u(A — 1) < 6m as m — > oo. (By this, we mean that if v m is the maximum possible 
value for given m, then limsup,^^ v m /m < 5.) 

We make the normalizing assumptions /(0) = 1, /i(0) = 1, g(Q) = 0, </(0) = 1. (We could also assume 
/' (0) = 0, but this is somewhat unnatural in the cases of interest.) 

In this context, we recall the following variant of the Biirmann-Lagrange theorem: 

Lemma 3.1. Let 4>{t) and tp(t) be formal power series, where ip(0) = 0, ip'(0) ^ 0. Then 

^{t) = Y,^m\ (3-2) 

0<i 

where 

«i = [t^m^m- 1 - (3.3) 

The advantage of this formulation for our purposes is that the dependence on j below is encoded in a single 
power series. 

Corollary 3.2. Let the coefficients aj and Ci be related by the formal power series identity 

5>^" = h{t)Y J c i f{t) m - i g{t)\ (3.4) 

j>0 i>0 

with f, g, h as above. Then 



c, 

3 



where 



Proof. We write 



and thus 



aij = [t°]*»'/»(t)- 1 /(*) i - m ff(t)- i 



= J2 a v a ji ( 3 - 5 ) 
tg'(t) tf'(ty 



(3-6) 



9(t) fit) 

IP =h{t)Y,<*i j f(t) m - i g{t)\ (3.7) 



i>0 



h{t)-'f{ty™V=^a i3 (^j ; (3.8) 
the formula follows immediately from the lemma. □ 
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In particular, 

Corollary 3.3. Let f , g, h be as above, and suppose 

j2*jt j = Ht) cifitr-'gity. (3.9) 

j 0<i<fn 

Then for all i > m, 

We will thus need an asymptotic analysis of the coefficients a>ij. In fact, until Section [| we will need only 
the cases i = m + 1 and i = m + 2. 

To make the asymptotic analysis tractable, we need some additional assumptions (summarized in the state- 
ment of Theorem |3.6| below). First, we assume / >: 0, / ^ 1; in all of our applications, / is derived from a 
weight enumerator or theta series, so this condition is automatic. Our second condition, that 1/g y 0, is less 
automatic, but is easily verified in all cases of interest. 



Define a function Lg(t) — tg'(t)/g(t). By Lemma 2.2 applied to the Laurent series 1/g, we find Lg'(t) < 
for positive t within the radius of convergence of 1/g. Thus if the equation Lg{t) = has a positive real root, 
it must be unique. We assume the root exists, and denote it by to- Note that if 1/g has radius of convergence 
r, and lim^,,- g(t) = 0, then l/g(t) is eventually increasing, and thus Lg(t) is eventually negative; since 
Lg(0) — 1 > 0, this implies that Lg{t) = has a positive real root. 

We finally assume that to is within the open disc of convergence of f(t), f(t)/h(t), and f'(t)/h(t), and that 
l/h(t) > on [Q,t ]. 

Lemma 3.4. Let f and g satisfy the above conditions. Then f(t)/g(t) has a unique local minimum on (0,io). 

Proof. Since f/g >: 0, we conclude that f(e s )/g(e s ) is strictly convex on (— oo, logto). Now, lim t ^ f(t)/g(t) = 
+oo, and thus f(e s )/g(e s ) is decreasing in a neighborhood of s = — oo. It thus remains only to show that 
f(t)/g{t) is increasing near to- 
We compute: 

t J t j$j = 9{t)-\f{t) Lg(t)f(t)). (3.11) 
At t , this becomes f'{t)/g(t) > 0. □ 

Denote this local minimum by t\. 
Lemma 3.5. Let f , g satisfy the above assumptions. Then there exists a unique point t' in (0,to) such that 

M) _ /(*o). 



g(t' ) g(to)' 

moreover, t <t\. 



(3.12) 
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Proof. Uniqueness follows from strict convexity, so it suffices to show existence. But f(t)/g(t) is continuous 
and decreasing on (0, t{\, and converges to oo at 0, so attains every value greater than f {t\) / g{t\) , in particular 
/(*o)M*o). □ 

The basic idea behind the proof of Theorem |3.6| below is that, to first order, the relation corresponding to 

Cm+2 7— T c m+1 (O.lcSJ 

9{to) 

is at j = and is positive for j > Lg(t' )m. Moreover, if we perturb the relation by subtracting a small multiple 
of c m +i, the relation becomes positive at j = 0, at the cost of slightly reducing the range of positivity. I.e., 
we can use Lemma [2.3] to obtain a contradiction for all relaxations of the desired bound when m is sufficiently 
large. 

Theorem 3.6. Let /, g, and h be convergent real power series satisfying the following hypotheses: 

(I) . f,h = l + 0(t), g = t + 0(t 2 ). 

(II) . f. l/g y 0. 

(Hi). Lg(t) := tg' (t) / g(t) has a positive real zero; let the smallest such zero be denoted t$. 
(iv). f, l/g, f/h, and f'/h have radius of convergence > to. 
(v). 1/h is positive on [0, to]- 
For each integer m > 0, let d m be the maximum of is (A — 1) where A ranges over power series >z of the form 

A(t) = h(t) c t f(t) m -*g(ty. (3.14) 

0<i<m 

Then limsup,,^^ d m /m < Lg(t' ), where < t' < to is such that f(t' )/g(t' ) = f(to)/g(to)- 

Proof. Choose < t 2 < t' ; we will show that for all sufficiently large m, d m /m < Lgfa). As this will hold for 
all choices of t 2 , the theorem will follow from the continuity of Lg. We also choose < £3 < t 2 . 

We consider the linear combination c m+2 — f(t3)/g(t3)c m +i, as a vanishing linear combination of the co- 
efficients of A(t). In particular, if we let ajfe) denote the coefficient of [P]A(t) in this linear combination, it 
will suffice to show that for all sufficiently large m, ctj(t^) > for j = and j > Lg(t 2 )m. Indeed, if a choice 
of A(t) existed with d m > Lg(t 2 )m, this would give a positive linear combination of nonnegative quantities (at 
least one of which is positive), equal to 0, a contradiction. 

We compute: 

aj(t 3 ) = a {m+2}j - jjp-atm+i)] (3-15) 
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where g(t) := t 1 g{t). We thus need to show that for all sufficiently large m, 

[t k ]F{t)G{t) m > (3.18) 

for k = m + 2 and < k < (1 — Lg{ti))m, where 

G(t)=g(t)~\ (3-19) 
FW^^f^^Uf^V (3-20) 

v; y ' \g(h) g(t)J dt\ g (t)J K > 

We note the following properties: 

(i) . F(t) and G(t) are power series with radius of convergence > to- 

(ii) . G(t) h 0. 

(iii) . F(t ) > 0. Indeed, the first two factors are clearly positive, while the second factor is positive since 

m = m < m (3 . 21) 

g(to) g(t' ) g(t 3 ) K ' 

Similarly, the third factor is positive, since to > t-y. 

(iv) . F(t) is positive on [0, t^). This time, the second and third factors are negative; we note the limit F(0) = 1. 



In other words, the hypotheses of Lemma 2.3 apply, taking / = [0, t2\ U [to, to + e] for appropriate e > 0. It 
follows that for all sufficiently large m, 

[P]F{t)G{t) m > (3.22) 

when j/m £ [0, 1 — Lg(t 2 )] U [1, 1 — Lg(t + e)]. For sufficiently large m, (m + 2)/m e [1,1 — Lg(t Q + e)], so we 

are done. □ 

Example. Let A(x,y) be the weight enumerator of a doubly-even binary self-dual code of length n. Then 
by Gleason's theorem, we have: 

A(x,y)= J2 c^ + UxV + y'r^ixVi^-y 4 ) 4 ) 1 . (3.23) 

0<i<[n/24] 

Defining a power series A(t) by A(t) — A(l,t 1 ' 4 '), we obtain: 

A(t)= c t (l + Ut + t 2 ) n / 8 - 3l {t{l-t) 4 ) 1 (3.24) 

0<i<[n/24] 

= (l + Wt + t 2 )"/ 8 - 3 ^ 24 ] J2 c l ((l + 14< + t 2 ) 3 )["/ 24 l- l (t(l-i) 4 ) 4 . (3.25) 

0<i<[n/24] 

We will apply the main theorem, with 

f(t) = (l + Ut + t 2 ) 3 (3.26) 

g{t)=t{l-tf (3.27) 

h{t) = (1 + Ut + £2)n/8-3[n/24] _ (3.28) 
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We clearly have f,l/g,h> 0, and each of /, 1/g, f/h, f'/h has radius of convergence at least 1; indeed, except 



for these are polynomials. Since limt— a g(t) = 0, the hypotheses of Theorem 3J3 are satisfied 



We easily compute to — 1/5; we then find that t' is the unique solution in (0, 1/5) of the quartic equation 

* 4 - 644t 3 + 6* 2 - 644*+l. (3.29) 

We thus obtain the bound 

limsup^i < In^L, (3.30) 
or, since v(A(t) - 1) = d(C)/4 and m = [n/24], 

d(C) l-5t' 

limsup-7-^ < — : (3.31) 



Ur** n(C) - 6(1 - *(,) 
In fact (as we will explain below), t' has the simple closed form 



1 - 5- 1 / 4 



4 



when substituted in, this simplifies (again explained below) to give 

limsupd„/n < i(l - 5~ 1/4 ). (3.33) 

4 Applications I 

We generalize the previous example as follows (compare the Gleason-Pierce theorem |HJ ) : 

Theorem 4.1. Let q and c be chosen with q > 1, such that either c = 1 or (q,c) <E {(2,2), (2,4), (3,3), (4,2)}. 
Let Ci be a sequence of formally self-dual codes over an alphabet of size q with all Hamming weights divisible 
by c; suppose furthermore that as i — ► 00, n(Ci) — ► 00. Then 

limsup ^ < - ( C + 1)-Vc). (4.1) 

i^oo n(Ci) q 

Proof. Given such a code C, of length n, let Ac(x,y) be its weight enumerator, and define a power series 
A(t) = A c (l,r^ c ). Then from the various Gleason theorems [|l8L Section 7], we conclude 

c=l:A(*) = (l + (V?-l)*) nmod2 Ci((l + (^/q- 1)^)^-^(1-0 (4.2) 

0<i<[n/2] 

(q,c) = (2,2) : A(t) = (1 + t )(n/2)mod4 ^ c .(( 1 + t )4)[n/8]-i( t ( 1 _ t )2y (4 3) 

0<i<[n/8] 

(g,c) - (2,4) :A(t) = (1 + 14f + * 2 )(™/ 8 ) mod3 ^ c,((l + 14* + t 2 ) 3 )[ n / 24 l- l (*(l - *) 4 )* (4.4) 

0<i<[n/24] 

(<Z,c) = (3,3) : A(*) = (l + 8*) (n/4)mod3 ^ c s ;((l + 8*) 3 ) [ " /12] ^(*(l - *) 3 ) 1 (4.5) 

0<j<[n/12] 

(q, C ) = (4,2) : A(t) = (l +3f )[n/2]mod3 Cl ((1 + 3*) 3 ) [n/6] ~ l (* (1 - i) 2 )\ (4.6) 

0<i<[n/6] 
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for appropriate coefficients Cj. In particular, we are in the scenario of Theorem 3.£, with g(t) = i(l — t) c . In 
each case, /, f /h and f /h are all clearly polynomials; since 1/g >z 0, has radius of convergence 1, and g(l) = 0, 
the hypotheses of the theorem are satisfied; it remains to compute Lg(t' ). 
Since 

Lg(t) . izil±£S, (4 . 7) 

we find to = 1/(1 + c). To compute £q from to> we proceed as follows. 

Define new series F(i) = f(t c ) and G(t) = g(t c ). From the MacWilliams identity and the fact that / and 
g are linear combinations of power series coming from weight enumerators, we find that there exists an integer 
no such that: 

F ( _ iz* ) = qW2 (1 + (g _ i) t )-»ojr( t ) (4.8) 



l + (g-l)t 



G 

Dividing these equations, we find 



( i + fa-nt ) = qn ° /2(1 + ( q - m-^Git). (4.9) 



In terms of / and g, this becomes: 



We thus conclude that 



f / / 1 - t x l c \\ f 

1 I > I » = % (4.11) 



3 I Vl + (9-l)« 1/c / / 5 



.l + fe-ljfc+l)- 1 / 

since we readily verify < <q < to- 

Similarly, to compute Lg(t' ), we differentiate the functional equation for G at 4jj , obtaining: 



l 9 (Q = - 4 /c ) + W f fa- 1 )^-^- 2 )-^^ = (i^lH (1 _ (4 . 13) 

cq y q J cq 

since by definition Lg(t ) = 0. Since m = [n/n ] and v{A{t) — 1) = <i(G)/c, we mutiply this by c/n to obtain 
the desired bound. □ 

Remark. In particular: 

(i) . For doubly-even self-dual binary codes, d/n < (1 - 5" 1 / 4 )/2 = .1656298476. 

(ii) . For self-dual ternary codes, d/n < (2 - 2 1 / 3 )/3 = .2466929834. 

(iii) . For even self-dual additive codes over GF(4), d/n < (3 - 3 1 / 2 )/4 = .3169872982. 

(iv) . For singly-even, formally self-dual binary codes, d/n < (1 - 3" 1 / 2 )/2 = .2113248655. 
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(v). For formally self-dual codes over GF(q), d/n < 1/2 — l/2q. 

Strictly speaking only the case (q,c) — (3,3) is new; the case (q,c) = (2,4) was shown (via a rather different 
proof) in 0), while in the remaining cases, the bound obtained is worse than the JPL || or Aaltonen §, § 
bound, as appropriate. But for singly-even self-dual binary codes, see Theorem 5.1 below. 



For even modular lattices (see [[T2|, 1 13 , |T{|; note that a unimodular lattice is 1-modular), we proceed 
similarly. We recall Dedekind's r\ function 

r)(z) = e" 2 / 12 Y[(l - e 2 ™) (4.14) 

i 

and the Eisenstein series 

E 2 (z) = -LAi og?7 ( z ) (4.15) 
Zm az 

= ^-E(E fc ) e2 ™" 2 - ( 4 - 16 ) 

l<m k\m 

Theorem 4.2. Let N be one of the integers {1, 2, 3, 5, 6, 7, 11, 14, 15, 23} ; and define 

E ( 2 N \z) = ^ mE 2 {mz). (4.17) 

m\N 

Then for any sequence Aj of even, strongly N -modular lattices of dimension tending to 00, 

v M(AQ Nz q 

hm sup — — -— — < — — , (4.18) 

i^oo dim(Ai) 2iti 

where z$ is the unique zero of E% on the positive imaginary axis. 

Proof. Let A be an even, strongly A^-modular lattice, with theta series 0a- Then (]l2|,|^3l) Q\{e^ lz ) can be 
written as a weighted-homogeneous polynomial in @A (e mz ) an d g{ e27Tlz ) — (Y[ m \N i]( mz )) 24 ^^~' mlN ™\ where 
Ao is the lowest-dimensional even A'-modular lattice. 

Clearly 9 Ao (i 1/2 ) h 0; using the product formula for 77, we also conclude that l/g(t) >z 0. Thus Theorem 
3.6 applies, and it remains only to compute Lg(t' ). 

We first compute, with t = e 27Tlz , 

Lg(z) = log J| )](mz) 24/E H» m ) (4.19) 

m\N 

24 -E { r\z), (4.20) 



Z~2m\N 

and thus t = e 2Trlz ° . 

We have the transformation laws: 



f(-l/Nz) = N n °/ 4 (z/i) n °/ 2 f(z) (4.21) 
g(-l/Nz) = N n °/ 4 (z/i) no/2 g(z), (4.22) 
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where 



no . "SOL. (4.23) 

We thus conclude that t' a = e - 2 ^/ Nz °, and that 

2:0^0 2 zoATn 

g{a) = + z ° L5(<o) = i^-- (4 - 24) 

Multiplying by 2 (since the lattices are even) and dividing by no (since m = [n/fto]) gives the required bound. □ 
Remark. Numerically, we have: 

N = 1 : -t^ttt £ .0833210664 N = 2 : - < .1246710056 

dim(A) dim(A) 

N = 3 : Jf^l s < .1643714543 AT = 5 :-t^ttt £ -2351529896 

dim(A) dim(A) 

N = 6 :-^4tt < .2414115212 AT = 7 : -^jj . < .2957105217 

dim(A) dim(A) 

N = 11 : r^rl \ £ -3973198712 AT = 14 :-^rr £ .4266498017 

dim(A) dim(A) 

AT = 15 :-^4tt < .3206725342 AT = 23 : M ^ < .6262824896 

dim(A) dim(A) 

Again, aside from N = 1, AT = 2, JV = 3, the obtained bound is worse than that implied by the Kabatiansky- 
Levenshtein bound on sphere packings 0. 

We finally consider self-dual codes over Z4. As in |Q, bounding the Hamming or Lee distance reduces 
to a consideration of the dual distance of doubly-even binary codes. At length a multiple of 8, the bound on 
self-dual doubly-even codes applies; for other lengths, we shorten the code up to 7 times, without affecting the 
asymptotic bound. We thus obtain: 

Theorem 4.3. Let C, be a sequence of self-dual codes over Z4, with length tending to 00. Then 

d H (Cj) 1 - 5- 1 / 4 

hmsup < (4.25) 

limsup^pl < 1-5- 1 / 4 , (4.26) 

i — >oo n\Lji J 

where dff(C) and dh{C) are the minimum Hamming and Lee weights of C , respectively. 

For the Euclidean distance, the situation is more complicated. We use the following lemma: 

Lemma 4.4. Let C be a Type II self-dual code over Z4 (all Euclidean norms divisible by 8), and let A(x,y,z) 
be its symmetrized weight enumerator. Define a power series 

A{t) := A(l, 2- 1/4 t 1/s {l + i) 1/4 , i 1/2 ). (4.27) 
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Then ^4^0, and 

v(A(t) - 1) = d E {C)/%. (4.28) 

Moreover, for appropriate coefficients C;, 

A(t)= c i (l + 60i+134t 2 + 60f 3 + i 4 ) (n/8) - 3l (i(l-<) 6 (l-i 2 ) 2 ) 1 . (4.29) 

0<i<[n/24] 

Proof. Consider a monomial x a y h z° of A(x, y, z); note that b + 4c must be a multiple of 8, and in particular b 
is a multiple nof 4. Under the specified substitution, this yields 

2-6/4(1 + t )6/4 t (6+4c)/8 j (430) 

a polynomial with nonnegative coefficients and with valuation (6 + 4c)/8; the first two claims are immediate. 

Now, the Gleason theorem for Type II codes over Z4 states that A(x, y, z) is a (weighted homogeneous) 
polynomial in 

(x + zf + (x - z) 8 



2 



128y 8 (4.31) 



i 6 = {x 2 z 2 (x 2 + z 2 ) 2 - Ay 8 ) ((a; 4 + 6x 2 z 2 + z 4 ) 2 - 64y 8 ) (4.32) 

24 

2 1 _2\ r>„.4\2 



y 8 {x 2 - z 2 f (4.33) 



h s = {xz{x z + z z ) -2 y y. (4.34) 
Under the substituion, we have: 

8 h-> 1 + 60i + 134i 2 + 60i 3 + t 4 (4.35) 

0i 6 i-> (4.36) 

024 ^ - t) 6 (l - t 2 ) 2 (4.37) 

/ig 1 ^ (4.38) 

The remaining claim follows. □ 

Remark. The above substitution is inspired by the proofs used in 0] and JT9[ ] , which involve lifting the code to 
a lattice and analyzing the resulting theta series. The resulting substitution takes the polynomial 

2y 4 - xz{x 2 + z 2 ) (4.39) 

to 0; solving gives the above substitution. 



We apply Theorem with f(t) = (1 + 60< + IMt 2 + 60t 3 + i 4 ) 3 , g(t) = t{l - t) 6 (l - t 2 ) 2 , to obtain: 
Theorem 4.5. Let Ci be a sequence of Type II self-dual codes over Z4 with length tending to 00. Then 

lim sup dB f. ^ < .3332625492 = 2 - 2x, (4.40) 
where x is the positive real root of the polynomial 

llx 16 + 2112a; 14 - 8525a; 12 + 15048a; 10 - 15218a; 8 + 9552a; 6 - 3718a; 4 + 828a; 2 - 81. (4.41) 
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Remark 1. The field Q[x] is that generated by 

2- 1 / 4 ^ /8 (l+t ) 1 /4 and ^/ 2 , (4.42) 

where t = 2v ^~ 3 . 

Remark 2. The Type II hypothesis is removed in Theorem |5.3| below. 

5 Applications II 



We now consider applications to which Theorem 3J3 does not directly apply, but for which the same basic idea 
can be used. 

We first extend the bound for doubly-even self-dual binary codes to general self-dual binary codes. The 



main idea is that, using Lemma 2.6, we can reduce positivity of coefficients of the form 

[V]Fi(t)F2(t) l G(1*) m (5.1) 
to positivity of Fi(t), for all sufficiently large / and all sufficiently large m; this because 

[t 23 ]Fi(t)F 2 {t) l G(t 2 ) m = [V] (F 1 (y/i)F 2 (Vt) 1 + F 1 (-Vt)F 2 (-Vt) 1 ) G(t) m (5.2) 

it^]F( t )F 2 ( t yG(t>r = m Gm (5 . 3) 

Theorem 5.1. Let C% be a sequence of self- dual binary codes of length tending to oo. Then 

d(d) 1-5- 1 / 4 , 

hmsup-^< . (5.4) 

Proof. We recall that a self-dual binary code C has associated to it two enumerators: its weight enumerator 
A(x, y) and its "shadow" enumerator S(x,y) ||, For appropriate coefficients Cj, we have: 

A( X ,y) = {x 2 +y 2 r™^ c^+y^t^-W^-y 2 ) 2 )' (5.5) 

0<i<[ra/8] 

S(x,y) = {2xy)^ od ^l 2 £ Cl {(2xy) 4 )^-\-{x 4 - y 4 ) 2 / A)* (5.6) 

0<i<[n/8] 

Here Cj can be characterized (up to a multiplicative constant) as the unique linear combination of the first i + 1 
nontrivial coefficients of A(x,y) that is also a linear combination of the first [n/8] — i nontrivial coefficients of 
S(x,y). 

We define power series A(t) and S(t) by 

A(t) = A(l,t^ 2 ) (5.7) 

= (l + t )(«niod8)/2 ^ Q ( ^ + ^4) [n/8] -i (j. ( j _ q2 y (53) 
0<i<[n/8] 

S(t) = t -(««w>d8)/8 S r^ lj t l/4j (59) 

= 2 (nmod8)/2 ^ Cl (16i)["/ 8 M(_(1 - i) a /4)\ (5.10) 
0<i<[n/8] 



1 Note that this includes the doubly-even case, for which S(x,y) = A(x,y). 
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If we let otij denote the coefficient of [P]A(t) in Cj, and let [3ij denote the coefhcient of [P]S(t) in Cj, Biirmann- 
Lagrange tells us: 

ay = [f- J ](l - 6t + i 2 )(l + ^""/a-i^ _ t j-M-i (5 n) 

= [t* _i ](l - 6* + i 2 )(l + f) 6i "" /2 (l - t 2 )- 24 - 1 (5.12) 
/% = (_i)*2 6l -™/ 2 [t["/ 8 ]-^](l + *)(l - i)-*- 1 . (5.13) 

Set to = 5~ 1//2 , = f ^5-1/4 ) ■ an d let to be an even integer of the form n/12 + 1/6. Now, consider the 
coefficient 

(l+*o) 4 K , A , 

Cm+1 - t (l-t y Cm - (5 ' 14) 

This is manifestly a negative linear combination of [P]S(t) for < j < [n/8] — to (so the same will be true for 
small perturbations). It thus suffices to consider the coefficients of [P]S(t). By the remark above and Lemma 
2.6 , this reduces to showing that 

is positive for < t < t' a , and can be perturbed to be positive at to as well; this is clearly the case. □ 

Similarly, the bound for even self-dual additive codes over GF{4) extends to the general case; the resulting 
bound is still weaker than the Aaltonen bound, however. A similar remark applies to formally self-dual binary 
codes, with shadow defined by S(x,y) = A(y,x); there the bound for singly-even f.s.d. binary codes extends. 



In that case, to ~ n/4, instead of the obvious analogue n/3; also we must apply Lemma 2.6 to the coefficients 
of S as well as to the coefficients of A. 
For lattices, we have: 

Theorem 5.2. Let N be one of the integers {1, 2, 3, 5, 6, 7, 11, 14, 15, 23}, and define 

E ( 2 N \z) = mE 2 {mz). (5.16) 

m\N 

Then for any sequence Aj of strongly N -modular lattices of dimension tending to oo, 

hmsup— — -— — < — — , (5.17) 

i^oo dim(Ai) 2m 

where zq is the unique zero of E% on the positive imaginary axis. 

Proof. As above; the case N odd is analogous to the self-dual binary code case, while the case N even is 
analogous to the formally self-dual binary code case. The only respect in which the proof is not straightforward 
(using the formulae of fl9| ) is in dealing with the "other" genera (not covered by Theorem 2 of In each 

case, direct summing by a suitable A-modular lattice places us into the "good" genera, and we can proceed 
from there; the only effect is to multiply the power series in question by a theta series, which clearly has no 
effect on positivity. □ 
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Our last shadow application is to codes over Z4: 
Theorem 5.3. Let Ci be a sequence of self-dual codes over Z4 with length tending to 00. Then 

limsup < .3332625492 = 2 - 2x, (5.18) 

where x is the positive real root of the polynomial 

llx 16 + 2112a; 14 - 8525a; 12 + 15048a; 10 - 15218a; 8 + 9552a; 6 - 3718a; 4 + 828a; 2 - 81. (5.19) 
Proof. We define 

A(t) = A(l, 2- 1 'H 1 / 4 {\ + t 2 ) 1 ' 4 , t) (5.20) 

S(t) = r ( nmod8 )/ 8 (l + t )-(nmod4)/4 S ( 1) 2 " 1 /4 t l/8( 1 + t )l/4 ) 

and observe that A and 5 are polynomials with nonnegative coefficients. Since 

A(t)= J2 c 4 ((l + i) 8 )"/ 8 - l (t(l+< 2 )(l-t) 4 r (5.22) 

0<i<[n/8] 

S(t) = (1 + *)[™/ 4 W°d2 ^ (-1) 1 (64<(1 + <) 2 )["/ 8 l- 4 ((l - t) 4 /8)\ (5.23) 

0<i<[n/8] 

we can proceed as in the case of self-dual binary codes. □ 

Our final application is to quantum codes. Formally, a g-ary quantum code corresponds to a self-orthogonal 
codes C over an alphabet of size q 2 ; the objective is to bound the minimum weight of the nonlinear code C 1 - — C. 
We recall the following (the nonbinary extension of (h], Theorem 6]): 

Lemma 5.4. Let Q be a q-ary quantum code of length n and dimension K . Then there exist polynomials 
C(x,y) and D(x,y), homogeneous of degree n, such that 

C { x,y) = c C + {q2 - l)V , X -^) (5.24) 

q q 

n „,x + (q 2 — l)y x — y, 
D(x, y) = -D( ^ y -), (5.25) 



and satisfying 



C(1,0) = 1, (5.26) 



v(C(l,t)-D(l,t)) =d (5.27) 
C(l,t)h^^(C(l,t)-D(l,t)) (5.28) 
C(l,t)-D(l,t)h0. (5.29) 
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In particular, if C and D satisfy the constraints for a given value K > 1, they satisfy them for all smaller 
K, including for K = 1. Thus if we replace ( 5.28 ) by the condition 

C(l,t)hO, (5.30) 

the resulting bound will apply to all quantum codes. 

Theorem 5.5. Let Qi be a sequence of q-ary quantum codes of length tending to oo. Then 

dm < i / _ i 

n{Qi) - 2 V <T 



l imsup ^™ < t (1-4) . (5.31) 



Proof. We consider the case of odd length n — 2m + 1; the case of even length is analogous. Setting C(t) :- 
C(l,t), D{t) = D(l,t), we observe that 



C(t) = (l + (q-l)t) J2 c i ((l + (q-l)t)*r- i (t(l-t)) i (5.32) 

0<i<m 

D(t) = (l-(q + l)t) ]T d id 1 + (? -l)*) 2 r"W -«))', (5.33) 

0<i<m 

for suitable coefficients Cj and c^. Let 7-y be the coefficient of [i J ]C(i) in (extending as usual to i > to), and 
let <5y be the coefficient of \V]D(t) in e?j. We find: 

7( m+ i), = [t m+2 - J ]{l ~{q+ l)t)t{\ - i)(l - t)-" 1 - 3 (5.34) 

7(m+2 W = [* m+2 - J '](l - (g + l)t)(l + (q- l)t) 2 (l - t)- m ~ 3 (5.35) 

<WDi = [i m+2 ^](l + (9 - !)*)*(! - - *) _m_3 (5-36) 

V+2)i = [i m+2_J '](l + (9 - 1)*)(1 + (9 - 1)*) 2 (1 - i)" m " 3 - (5-37) 

If we instead expand Cj and dj in terms of [P]C(t) and [i J ](C(t) — D(t)), we obtain coefficients: 

(coeff. of [t j ]C(t) in c,) = 7ij (5.38) 

(coeff. of [V]{C{t) - D{t)) in a) = (5.39) 

(cocff. of [t°]C{t) in di) = 5ij (5.40) 

(coeff. of [t J ](C(t) - D(t)) in d;) = -S tj (5.41) 

We need a relation that is a nonnegative linear combination of the coefficients [P]C(t), positive at j = 0, as 
well as a nonnegative linear combination of the coefficients [P]C(t) — D(t) for j larger than the bound. 
Now, consider the relation 

(q + l)( Cl - (q + l) 2 c ) - (q- l)(di - (q + lfd ). (5.42) 
On [t j ]C{t) - D(t), this has coefficient: 

(q - l)r +2 - J ](l + (<? - 1)0(1 - 2t)(l - (9 2 + l)t)(l - i)-'^ 3 , (5.43) 



21 



which to first order is positive for j/m > 1 — 4f. Similarly, on [P]C(t), this has coefficient: 

2 [t»"+a-i](i - 2t)(l - ( g 2 + 1)<) 2 (1 - f)- m - 3 (5.44) 

which is positive for < j < m + 2 except in a neighborhood of j = and j = (1 — ■^)m. The construction of 
a positive perturbation is straightforward. □ 

Similarly, taking the shadow constraints Jl7[ into account, we obtain: 
Theorem 5.6. Let Qi be a sequence of binary quantum codes of length tending to oo. Then 

d(Q,) 3->/3 

limsup — — < — . (5.45) 

i-^oo "(Qi) 4 

Note that although as we have remarked, this is slightly worse than the Aaltonen bound, this is still a new 
result; in the quantum case, the Aaltonen bound is only known for a set of rates bounded above g], 

6 Extensions 

As we have remarked, many of our bounds are weaker than the appropriate "universal" bounds (JPL jnj; 
Aaltonen Q , j| ; Kabitiansky-Levenshtein || ; see Q for a survey) that hold even for non-self-dual codes of rate 
1/2 and non-lattice packings of appropriate density. Since others of our bounds are quite a bit stronger than 
the corresponding universal bounds, this strongly suggests that in no case is either bound tight for self-dual 
codes. The question then becomes that of how to improve the above bounds. 

We restrict our attention to the situation of Section |^; we will comment on the shadow and quantum cases 
at the end. 



Thus, let /, g, h satisfy the hypotheses of Theorem BJ: . As the above bounds resulted from considering the 



two relations c m +i, Cm+2 m place of the single relation c m +i, the obvious thing to try is a linear combination of 
c m +i for 1 < i < /. Somewhat surprisingly, this does not give rise to any improvement in terms of limsupd/m. 
We can see this as follows. The coefficient of [t J ]A(t) in such a linear combination will have the form 



[t l Wt)) P[ g (t) }t dt{g(t) 



~g(tr m , (6-1) 



(6.2) 



for some polynomial p. The point, then, is that since 

dt \g(t)J 

has opposite signs at t = to and at t = tg, while 

{m) (8 - 3) 

both have the same sign at the two points (the same value in the latter case), the corresponding coefficients 
will, to first order, also have opposite sign. In particular, we will never obtain a bound on d/m better than 
Lg{t'o). (This tends to explain why Krasikov and Litsyn obtained the same bound in the doubly-even binary 
case, despite a rather different argument, and the unlikelihood that the bound is optimal.) 
On the other hand, we do obtain a slight lower-order improvement: 
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Theorem 6.1. Let f, g, h, d m , to, t be as in the hypotheses of Theorem 3.6, and suppose further that 



is smaller at to than at t . Then 



lim sup ■ 



\dtW) ) 
-tLg'(t) 

■, - L 9(t'o) m 



(6.4) 



(6.5) 



Proof. For each k > 0, let be the relation corresponding to the polynomial p(t) = (t — /(to)/ g(to)) k above; 
that is: 



„(fc) 



EM) 



0<i<k 



/fa) 

g(t ) 



(6.6) 



Thus the coefficient of [t j ]A(t) in c (fe ) is 



Also, define coefficients by 



" 3] (M*))'*(s(*)) ( ff(*) sfo)^" 



(6.7) 



,(fe) .. 



[t™^].g(t)- 



We will consider relations of the form 



J2 b k m k ' 2 e^ 



(6.8) 



(6.9) 



0<fc<2 



with bi > 0. We first claim that for any such relation, the coefficients of [t 3 ]A(t) are nonnegative when 
j > Lg(t' )n + n fc /( fc+1 ) and n is sufficiently large. Indeed, in a neighborhood of j/n = 1, this follows from the 
estimate of Lemma 2.3; the terms for k < I are o(m ! / 2 ), while the term for k = I is fl(m 1 ^ 2 ) and positive. In 
the other region, we use the 0(n _1 ) estimate of Lemma 2.4, in which the main term is positive and of order 
jl( n - fc /(fc+i)) > o( n -i). 



In the remaining neighborhood of t' , we use Corollary 2.5. We thus have: 



a L g (t> )m+x ^ = C[(C' 2 ) k+1 (C' 3 m)- k / 2 h k (x/^a 3 ) + o(m- k ^ max(|x| fc , 1)), 



where 



c{ = MO -1 > 

c - t > - (M\ 

C 3 = -t' Lg'{t' ) > 0. 



> 



(6.10) 

(6.11) 
(6.12) 

(6.13) 
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In particular, if x is bounded above the largest zero of 

J2 b k C[{C' 2 ) k+ \C',)- k l 2 h k {x/^Ci), (6.14) 

0<fc<( 

then we have positivity for j > Lg(t' )m + x^/m; when j < Lg(t' a )m + m~ l ^ l+1 \ the error term is uniformly of 
smaller order than the main term. 
Similarly, 

4 fe) = C 1 (C 2 ) fc+1 (C 3 m)- fe / 2 ^(0) + o(m- fe / 2 ), (6.15) 

with 

d = hito)- 1 > (6.16) 

C 3 = -t Lg'{t ) > 0. (6.18) 

Thus given any choice of b k such that 

q:= b k C\(C 2 ) k+1 (C 3 )- k ^h k (0)>0, (6.19) 

0<fe<i 

if x(b) is the largest zero of 

p(x) := ^C[{C' 2 ) k+ \C' 3 )- k ^h k {x/^/Ci), (6.20) 
o<fc<; 

we have the asymptotic bound 

d m - Lg(t' )m 

lmisup — — ^ x \b)- (6-21) 

m^oo TTl ' 

To construct a good relation, we will need some further properties of the Hermite polynomials, all classical 
results: 

Lemma 6.2. The polynomials h k (x) are the unique monic polynomials such that 

-= hj{x)h k {x)e- x2 l 2 dx = 5 jk k\. (6.22) 

V27T J-oo 

Furthermore, we have the three-variable generating function: 

, „ , 1 / t 2 x 2 /2-txy + t 2 y 2 /2\ 

£ h k {x)h k {y)- = -j== exp ^ j , (6.23) 

convergent for \t\ < 1. Finally (Christoffel-Darboux): 

h k (x)h k -i(y) - h k _i(x)h k (y) = x - h 1 (x)h j (y) 

x-y _ ^ ?! ^ ' ' 
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Now, given p(x), we can compute b k using orthogonality; we find: 



C[C 2 



i i 



J p{xyfa.)h k {x)e- x2 l 2 dx] 



using the same formula to define b k for k > I gives b k = 0. But then q can be computed as 

C X C 2 1 



C[C 2 y/2ff 



where 



J pix^C^y^)^' 2 ^, 



C^M\ h k (x)h k (0) 



k\ 



Using the three-variable generating function, we find: 



e- x2 / 2 y{x) = e-* 2 ' 2 ]T h k (x)h k (y)- = CT 1 cxp(-(x/C) 2 /2), 



fc>0 



where 



^(c>vc- 3 ) 2 -(c 2 ^c>y 



C :-- _ 
and we must satisfy the additional requirement 



precisely the additional hypothesis above. But then 

_ C1C2 1 
C[C 2 \p2/K 



J P (xC^C r 3 )e- x2/2 dx. 



Since 



C\C<2, 

C[C 2 



<o, 



(6.25) 



(6.26) 



(6.27) 



(6.28) 



(6.29) 



(6.30) 



(6.31) 



(6.32) 



the problem reduces to the following. For an integer k > 0, how small can we make the largest zero of a 
polynomial p{x) of degree 2k + 1 subject to the condition 



J p(x) 



e~ x2 l 2 dx < 0? 



Let x^ be the smallest zero of h k (x), and consider the polynomials 

Px(x) = h k (x) 
h k (x) 



P2(x) = 



(fe) 



(6.33) 

(6.34) 
(6.35) 
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By Christoffel-Darboux, we compute 



M X ) = h k . 1 (4 k) r 1 E hAx)h ^ )] (6.36) 

Q<j<k ■'' 



We thus have the following integrals 
1 



OO 

Pl {x)p2(x)e~ x2/2 dx = (6.37) 



i r°° h ( wy^ 

-L / p 2 (x) 2 e-* 2 / 2 dx = h^ix^)- 2 E ^r 1 - > 0. 
V2W-OC < J<fc J- 

Now, consider the polynomial 



(6.38) 



p(x) = Pl (x)p 2 (x) - ep 2 {xf + e 2 = (x - xf ] - e)p 2 {x) 2 + e 2 (6.39) 

for e > small. This certainly satisfies the negative integral condition; on the other hand, its only zero is that 

(fe) 

near Xq . 

Thus when I — 2k, we obtain a bound of the form: 

lim sup d ™ L ^) m < cyc^f . (6.40) 

m — >oo V W 

Since 

x^ fe) = -2\Zfc + 0(A:- 1/6 ), (6.41) 

we are done. □ 

Remark 1. The additional assumption is satisfied in all of the applications above; this is a trivial calculation 
for all but the lattice cases, in which it follows from the transformation law. 

Remark 2. The fact that taking 2k coefficients gives an improvement proportional to 2\J km suggests that to 
obtain a first-order improvement, we will need to let k grow linearly with m. 

Remark 3. The involvement of the polynomial h^{x) 2 /{x — Xq ) is somewhat reminiscent of the approach of 
Levenshtcin j^] to the universa bounds, in which different orthogonal polynomials occur. There the polynomial 
giving the bound is 

{Pk(x)p k -i(y) - p k (y)Pk-i(x)) 2 /(x - y), (6.42) 

which reduces to our polynomial when pk — and y — x^ . While in those bounds, one lets k grow linearly 
with n, one can also consider finite k; in that case, the improvement is f2(y / n), as in our case. 

For the shadow cases, roughly the same argument applies; for instance, in the self-dual binary code case, we 
take relations starting with m = n/12 + (logn) 2 /6. Positivity near j = m is immediate (since multiplying by a 
positive power series leaves positive initial coefficients positive) ; in the remaining region, we find that restricting 
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the integral to an interval \x\ < n~ e and replacing (1 + i)( logn ' 2 by (1 + |i|)( log ™) 2 gives negligible relative error. 
The argument then proceeds as before. 

For the quantum cases, the difficulty is in choosing the relation. Basically, one defines and as above, 
and considers a linear combination of 

m k ' 2 cf and m^' 2 {{q + l)c { k) - (q - l)d ( k) ). (6.43) 

The first set of relations has no effect on the coefficients of [t J ']C(i) — D(t), so as above, we essentially obtain 
an arbitrary polynomial here. Similarly, they have a lower-order effect on the coefficient of [t°]C(t); we thus 
end up with the same constraints on this polynomial as above. On the other hand, near j/m = Lg(t' ), the 
relations have the same order behavior; in this neighborhood, we may thus choose an arbitrary polynomial of 
degree k + 1, so have no additional constraints. 
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